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Abstract. The Mallows model on Sn is a probability distribution on permutations, q'^^^''') /Pn{q), where d(-K, e) 
is the distance between tt and the identity element, relative to the Coxeter generators. Equivalently, it is the 
number of inversions: pairs where 1 < « < J < n, but n-i > ttj. Analyzing the normalization Pn{q), Diaconis 
and Ram calculated the mean and variance of d(7r, e) in the Mallows model, which suggests the appropriate 
n — > oo limit has q„ scaling as 1 — f3/n. We calculate the distribution of the empirical measure in this limit, 
u{x, y) dx dy = lim„^oo ^ 5Zr=l ^(i,Tri)- Treating it as a mean-field problem, analogous to the Curie- Weiss model, 
the self-consistent mean-field equations are g^g^ In u{x, y) = 2/3u(x, y), which is an integrable PDE, known as the 
hyperbolic Liouville equation. The explicit solution also gives a new proof of formulas for the blocking measures 
in the weakly asymmetric exclusion process, and the ground state of the Wq(0[2)-symmetric XXZ ferromagnet. 
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1. Introduction and Main Results 

The Coxeter generators of the symmetric group Sn are the transpositions (1, 2), (2, 3), . . . , (n — 1, n). The 
height of a permutation is defined distance to the identity element e, 

d(7r, e) = minj/j > : 3 ri, . . . , Tfc S {(1, 2), . . . , (n — 1, n)} such that tt = ti • • • Tfe} . 

More generaUy, d(7ri, 112) — d{n2^Tri, e). It is easy to see that c?(7r, e) = d{TT~^ , e). In fact, another formula is 

d(7r, e) = G : I < i < j < n and > nj}\ . 

The pairs («, j) are called inversions of tt. In [5], Diaconis and Ram studied the Mallows measure, which is a 
probability measure on Sn given by 

d(7r,e) 



Pn{q) ' 

with Pn{q) being a normalization constant. Actually, Diaconis and Ram studied a Markov chain on Sn for which 
the Mallows model gives the limiting distribution. This was followed up by another paper on a related topic 
by Benjamini, Berger, Hoffman and Mossel (BBHM) [3] who related the biased shuffle and the Mallows model 
to the asymmetric exclusion process and the "blocking" measures (of Liggett, see [TT], Chapter VIII, especially 
Example 2.8 and the end of Section 3). They did this using Wilson's height functions [20]. We will discuss 
this more in Section |HJ For now, let it suffice that Diaconis and Ram identified the explicit formula for the 
normalization which they remarked is the "Poincare polynomial" : 

- 1 



Further references for the statistical applications^ of the Mallows model can be found in their paper. 



Date: March 2, 2009. 

^Independently, a similar g-deformed combinatorial formula was explained for a problem in quantum statistical mechanics, the 
ground state of the ferromagnetic Wq(B[2)-symmetric XXZ quantum spin chain, by Bolina, Contucci and Nachtergaele [4]. We will 
comment more on this in Section 8. 
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Note that, physically speaking, one would define a Hamiltonian energy function i7„ : S'„ -^^ R as 

^''^^^ = ~J 1(0,00) (tt^ -TTj). 

l<i<j<n 

In this case, one thinks of vr = (tti , . . . , 7r„) as some type of constrained spin system, where each of the components 
TTi, . . . , 7r„ are spins in {1, . . . , 71} as in a Potts model. The choice of the normalization of the Hamiltonian is 
then standard for mean- field models. One would be most interested in the free energy 



UP) = In E - 



/3i/„(7r) 

I3n 



7reS„ 

For our purposes, we prefer to consider the mathematically simpler "pressure" 



Pn{f3) = - In (1 E 



g-,3i/„(7r) 



(Note that, contrary to the usual conventions of statistical physics, we have divided the partition function, which 
is Zn{(3) — X^TTinS e~^^"^'^\ by the infinite-temperature partition function Z„(0) = n!, which is equivalent to 
starting with a normalized a priori measure rather than counting measure on S'„.) It is trivial to see that this 
is given precisely by the Poincare polynomial described by Diaconis and Ram: 

1 P„(e-^/("-^)) 1 
Pn[p) = - In j = - In j . 

With this scaling, it is also easy to calculate the limit: 



p(/3) = lim p„(/3) = / In dx , 

n^oo _/p \ f3x J 

(which can be solved explicitly using the polylogarithm function). From this one can calculate the mean and 
variance in the limiting Gibbs measure. For instance, one can calculate E{d{TT, e) — n{n — l)/4)^ ^ v? 112 in the 
uniform measure on Sn- 

To go beyond the statistics of (i(7r, e) it seems worthwhile to study the empirical measure of tt: 

1 " 
n ^ — ' 

which is a normalized measure on {1, . . . , n} x {1, . . . , n}. More specifically, this is a random measure. Rescaling 
the discrete cube {1, . . . , n} x {1, . . . , n} to [0, 1] x [0, 1], it is easy to see that the random empirical measure 
converges, in probability, to the non-random Lebesgue measure, when /3 = (the uniform case). Our main 
theorem generalizes this result. 



Theorem 1.1. For any /3 e K, 

, f 1 /• 1 

>e = 0, 



lim lim F,\-'3/" J iy^ j(i/„^7r,/n) - / 



f{x,y)u{x,y)dxdy 

[0,1] X [0,1] 



for every continuous function f : [0,1] x [0, 1] ^ IR, where 

(/3/2)sinh(/3/2) 



i{x,y) = 



(6/3/4 cosh(/3[a; - y]/2) - e-PI^ cosh(/?[a; + y - l]/2))^ 



Note that (one can show) the limit — > gives 1. The proof of Theorem 11.11 uses a rigorous version of 
mean-field theory, as in the solution of the Curie- Weiss model. An interesting feature is that the self-consistent 
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mean-field equation leads us to the characterization of u as the solution of an integrable PDE 

^--^\nu{x,y) = 2f3u{x,y). 

It is not unusual for mean-field problems to lead to integrable PDE's. We demonstrate this briefly in the next 
section with the ubiquitous toy model, the Curie- Weiss ferromagnet. 

2. Toy Model: The Curie- Weiss Ferromagnet 

We include this section merely to point out that mean-field problems often do lead to integrable PDE. However 
the issue is serious: in fact there is a recent paper by Genovese and Barra which we recommend for more details 
|13j . Our approach merely summarizes their results (in our own words) as well as the earlier paper by Barra, 
himself 2J. The configuration space of the CW model is = {+1, —1}^ = {a = (cti, . . . , cr„) : cti, . . . , cr„ — 
±1}. For technical reasons, we choose the Hamiltonian as 

t ^ Af 



we assume t > and x G 



Therefore, defining 



we easily see that 



and 



where 



27V 
r-i Y-^ 



Defining mjq(a^ = N ^ X]i=i '^i^ which takes values in [—1, 1], we see that 
Hj^ ^ -n(*-^+ xniN 



Pn 



V 2 



M{(7,t,x) 



1 



-Q^PN{t,x) = N {{m%) - {ruNf) , 



(/) — {f)N,t,x — 



E 



/We- 



E 



d_ 

dx 



PN{t,x) , 



Actually it is easier to consider the "order parameter," 

UNit,x) = (mN)N,t,x 

from which pN{t, x) can be calculated by solving the ODE: 

for X e M, 
as a; — !■ ±oo. 

Then we see that uiq{t, x) satisfies the viscous Burgers equation (with velocity equal to the negative amplitude): 
^ UN{t, x) = UN{t, x) ^ UN{t, x) + 2^ • ^u^it, x) for t > and a; G M, 



■^PN{t,x) = UN{t,x) 
PN{t,x) - \x\ -> 



Uff{0,x) = tanh(x) 



for x G 



This is an integrable PDE, using the Cole-Hopf transform. See, for instance. Chapter 4 of Whitham, [TO] . 
Actually, this leads to a solution in terms of Gaussian integrals. The analogous transform in spin-configuration 
notation is the Hubbard-Stratonovich transform: 
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which "hnearizes" the dependence of the Hamiltonian on mAr, in the exponentiaL This trick is used to solve the 
Curie- Weiss modeL See, for example, Thompson |T5]. 

Note that in the N oo limit, one obtains u{t, x) = limAr^oo UN(t, x) being the vanishing-viscosity solution 
of the inviscid Burgers equation. Shocks correspond to phase transitions. The Lax-Oleinik variational formula 
for solutions of hyperbolic conservation laws applies. See for example. Section 3.4.2 of Evans [11]. In this context 
we claim that this is equivalent to the Gibbs variational formula, in the mean- field limit. We review this next. 

3. Gibbs Variational Formula 

Let us begin by considering a general problem in classical statistical mechanics. Suppose that A" is a compact 
metric space, and suppose that there is a two-body interaction 

h:X X X ^ MU {+CX)} . 

We assume that h is bounded below. Then for each N > 2, one can consider the mean- field Hamiltonian 
Hn ■■ ^ MU {-Hoo} 

Hn{xi, . . .,Xn) ^ jy _ X! h{x^,Xj). 

l<i<j<N 

Suppose that there is also an a priori measure /io on A", which we assume is normalized so that 

dna{x) = 1 . 

Then the thermodynamic quantities arc the partition function, 

ZnW) = f e-''^"(----")dA.o(a^i)---dA'o(a^iv), 



the pressure. 



and the Boltzmann-Gibbs measure 



^-0HNixi,...,XN) 

dii'li{xi,...,XN) = Zn{I3) diXo{xi)---d^io{xN). 

Physically, it is more correct to consider the free energy rather than the pressure, fN{P) — —^PNiP)- But we 
will consider pn{(3), which seems slightly easier to handle, mathematically. 

We will write ji^ for the measure d^^ (xi, . . . , xat) = (i/io(a^i) • ■ ■ dfio{xN)- Also, if / is a function, then we 
use the short-hand for J fd^. Then, according to the Gibbs variational principle, we have 

(1) pj^{p) = max ^[SN{fJ.N,fJ.f'^) - f3fJ,N{HN)] , 

Miv6A< + ,i(A'™) 1\ 

where SNifJ-N , fJ-f^) is the relative entropy (and A^+^i(A'^) denotes all Borel probability measures on X^) 

IJ,Q {(t){diJLM I d^^)) \i is absolutely continuous with respect to /i^, 
-oo otherwise. 



and (f>{x) = — a;ln(a;), which is if a; = 0. Also, the unique /ijv maximizing the Gibbs variational formula (the 
"arg-max") is the Boltzmann-Gibbs measure /i^. 

A natural ansatz for the optimizing measure is = , for some measure fi S A4+^i{X). Probabilistically, 
this means that all the xi, . . . , xn are independent and identically distributed. Technically, this cannot usually 
be exact for finite A^. But it leads to a simpler formula because 

SN{pi^,fio) = NSi{pi,fio) and fi^iH^) = ^^'(/i), 
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and one hopes that the formula may become exact in the thermodynamic Umit. Mark Fannes, Herbert Spohn 
and Andre Verbeure proved that this approach is rigorous in the N ^ oo hmit [T^ : 

Proposition 3.1 (Fannes, Spohn, Verbeure 1978). The limiting pressure exists, p{(3) ~ limpf^ooPNiP), and 
solves the variational problem 

p{f3) = max [5i(^,/io) - -/^^(/i)] . 

Moreover, any subsequential limit of the sequence (/i^) is a mixture of infinite product measures for /i's 
maximizing the right-hand- side of the formula above. 

Remark 3.2. Note that the Gibbs variational principle fip is true in general for all Hamiltonians whether they 
are mean-field or not. (See, for instance, Lemma II. 3.1 from Israel's monograph [15j . or any other textbook on 
mathematical statistical mechanics, for a rigorous proof which also applies directly in the thermodynamic limit.) 
But the product ansatz which seems to yield the formula from the proposition is not generally valid, since there 
are nontrivial correlations in the true Boltzmann- Gibbs state. Nevertheless Fannes, Spohn and Verbeure proved 
the mean-field limit in the N ^ oo limit, using de Finetti 's theorem ( which states that all infinitely exchangeable 
measures are mixtures of product states ) and properties of the relative entropy. 

Because one has — fix p,, one replaces the linear form pNiH^) by the nonlinear one fJ,'^{h). Fannes, Spohn 
and Verbeure actually proved their theorem more generally for quantum statistical mechanics models, such as 
the Dicke maser, but it also applies to classical models. For the quantum models, one replaces de Finetti's 
theorem by the non-commutative analogue, Stormer's theorem. (See pQ and references therein for a detailed 
survey of de Finetti's theorem, and refer to Fannes, Spohn and Verbeure's paper and references therein for the 
noncommutative analogue, which we will not need.) With Eugene Kritchevski, we tried to find a simpler proof 
of the specialization of Proposition 13. II to the classical case. But there were several errors in our proof, which 
have been brought to my attention by Alex Opaku, to whom I am grateful. Fortunately, Fannes, Spohn and 
Verbeure's original paper definitely does also apply to classical models. 



4. Application to the Mallows model 

We take for X the unit square [0, 1] x [0, 1]. Suppose that f,g : [0, 1] ^ R are probability densities: f,g>0 
and Jp fix) dx — g{y) dy = 1. For simplicity, later on, we also assume that there are constants < c < C < oo 
such that c < f,g < C. Then we take the a priori measure to be 

dno{x,y) = f{x)g{y)dxdy . 

We take the interaction to be 

h{{xi,yi), {x2,y2)) = d{xi - X2}0{y2 - yi) + 9{x2 - xi)6{yi - ^2) , 
where 6* : M ^ R is the Heaviside function, 



0{x) 



if a; > 0, 
if a; < 0. 



Since fiQ is absolutely continuous with respect to Lebesgue measure, all xi, . . . , xn and yi, . . . , yjy are distinct, 
with probability 1. (This is why we do not bother to specify 9 at the discontinuity point 0.) 

Let Xi < • • • < Xn and Yi < • • • < Yn be any points. Then for any cr, r e Sn, the symmetric group, we have 



dfiS 
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where is the Mallows measure on Sjy. So studying the limit of the /i^'s gives us direct information on the 
limit of p]y _ Por any fixed a G Sn, the permutation a^^r is uniform on Sn, if t is. Because of this, we 
have the following result for the marginal of /i^ on (xi, . . . ,xn), 

U{xi, . . . ,XN)dfi^{{xi,yi), . . . ,{xN,yN)) = / U{xi, . . . ,XN)f{xi) ■ ■ ■ f{xN)dxi ■ ■ ■ dxN , 

X" "'[0,1]" 

and the marginal on (yi, . . . ,yN), 

Uiyi,...,yN)dfiff{{xi,yi),...,{xN,yN)} = / U{yi, . . . ,yN)g{yi) ■ ■ ■ g{yN) dyi ■ ■ ■ dy^ , 

"'[0,1]" 

for all bounded, continuous functions U : (0, 1)^ R. Enforcing these conditions on the marginals, Proposition 
13.11 vields the following: 

(2) p{(3) = max Mo) - /^m' , 

MeA^+,i(/,3) 

where is the set of all probability measures ^ G such that d^{x,y) has marginals f{x)dx 

and g{y)dy. 

Suppose that is any arg-max of the right-hand-side of Since we have chosen to be absolutely 
continuous with respect to Lebesgue measure on X ^ the same must be true of /i. Otherwise the relative entropy 
would be — cxD. So we can write 

d^[x,y) — u{x,y)dxdy. 
Then it is easy to see that the Euler-Lagrange equations for ^ are 

(3) \TLu{x,y) = \tl f{x) +\TLg{y) + C - (3 / u{x' ,y')[e{x ~ x)e{y' - y) + e{x' - x)0{y - y')]dx' dy' , 



IX 

for some constant, C < oo. Therefore, u solves the equation 

u{x,y) = ± f{x)g{y)e~'^ J'^ ''^(''■y^'^''' 'y"»^'^' ^y"> '^y' for {x,y) S X, 

(4) { uix, y) dy = f{x) for x e [0, 1] , 

/o^ u{x, y) dx = g{y) for y g [0, 1], 

where is a normalization constant. 

Since u{x, y) solves an integral equation it can be differentiated both with respect to x and y. Doing so yields 
the partial differential equation 

(5) 

dxdy^^^^''^^^^ 

known as the hyperbolic Liouville equation. This equation arises naturally in differential geometry, related to the 
problem of choosing a metric on a given manifold. I am very grateful to S.G. Rajeev for important information 
regarding this PDE. One of the facts he imparted is the symmetry of the differential equation under the following 
general transformation: 

v{x,y) = F'{x)G'{y)u{F{x),G{y)) 

lnu(x,?/) = -——u[F{x),G{y)) 



dxdy dxdy dy \ F'{x) J dx \ G'{y) 

= 2pF\x)G'(y)u{F{x),G{y)) 
= 2/3i;(x,y). 
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So, if ^Ig^lnw = Pu then the same is true for v{x,y) = F'{x)G'{y)u{F{x),G{y)). 

Our real goal is to solve the Euler-Lagrange equation ([4]). But as a first step, we want to consider the Cauchy 
problem for ([5]). In other words, we want to consider the problem 



(7) 



gfg^lnu(a;,y) = 2l3u{x,y) for {x,y) G [0, Li] x [0,L2], 
u{x,0) = (f){x) for X £ [0, Li], 

u{0,y) = for y G [0,2.2], 



for some Li, L2 > and : [0, Li] ^ R, : [0, L2] M both positive and continuous. 

Note that gfg^ is a wave operator, with characteristics directed along x and y. Specifically, defining ^ = 
(x + y)/ \/2 and C, = [x — y)/ V2, we have gfg^ = i (^'^ ~ ^^)' ^^"^ usual wave operator. Therefore, D'Alembert's 
formula for solutions of the wave equation allow us to reformulate ([7]) as an integral equation, 

(8) \nu{x,y) ln(/;>(a;) + lni/'(y) - Ina + 2/3 / u{x\y') dx' dy\ 

J[0,i;]x[0.y] 

which we prefer. This equation is supposed to be solved for all (x, y) E [0, Li] x [0, L2]. We have introduced the 
number a = (/)(0), which we also assumed equals ip{0), for consistency since both are supposed to give u(0, 0). 
(Note that the initial surface, ([0,Li] x {0}) U ({0} x [0,^2]), is not a non-characteristic surface. This is the 
reason that our Cauchy problem does not require initial data for the tangential derivative of u even though 
the wave equation is second order.) We refer to Evans textbook for PDE's, (especially Section 2.4 on the wave 
equation and Section 4.6 on the Cauchy-Kovalevskaya theorem). 

As we will see, the symmetry ([6)) is the key to solving both the Euler-Lagrange equation ([4]) and the Cauchy 
problem ([5]). 

5. The Cauchy Problem 
We start with uniqueness for the Cauchy problem. 

Lemma 5.1. For any Li,L2 > 0, the Cauchy problem (0) having (j) ^ tp ~ a ~ 1 has at most one solution in 
the class of nonnegative integrable functions. 

Proof. Since cj) = ip — a ^ I, equation ([8]) simplifies to 

\nu{x,y) = 2/3 / u{x' ,y') dx' dy' . 

JlO,x]xlO,y] 

Assuming that u is nonnegative and integrable, this implies that Inn is bounded and continuous. Then, using 
these properties in the right-hand-side of the equation again (similarly as one does to prove elliptic regularity) 
we deduce that Inu is continuously differentiable and globally Lipschitz. In particular, it is continuous up to the 
boundary. 

Now suppose that there are two solutions u and v. Letting z = Inti — Inu, we have 

z{x, y) = 2/3 /" /"V - e-'-'^'^'^y'^Hx', y') dx' dy' . 
Jo Jo 

Since both Inu and Inv are bounded, we see that z is as well. Therefore, there exists a constant K < 00 such 
that |1 — e^^l < K\z\ for all values of z in the range. So we have 

\zix,y)\ < pK\\u\\^ f r \z{x',y')\dx' dy' . 
Ja Jo 

A version of Gronwall's lemma then implies that z = 0. We outline this now, although our argument can 
probably be improved. 
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Let Z{t) = svLp{\z{x,y)\ : {x,y) e (0,Li) x (0, L2) , xy < t}. Then we obtain, after making the change of 
variables {x,y) >—>■ {x,t) where t — xy, and using Fubini-ToneUi to integrate over x first, 

Z{t) < PK\\u\\oo [ \n{t/t')Z{t')dt' . 
Jo 

We rewrite this as ^ 

Z{t) < l3K\\u\\oo I Mt/LiL-2)-Ht'/LiL2)]Z{t')dt' . 
Jo 

Since \n{t/ L1L2) < for t < L1L2, and since Z > 0, we can drop the term ln(t/LiL2)Z{t') in the integrand to 
obtain ^ 

Z{t) < l3K\\u\\^ f \\ii{t'/LiL2)\Z{t')dt' . 

JO 

Finally, setting ({t) = \ \n{t' / LiL2)\Z{t') dt' , this leads to 

C'(i) < l3K\\u\\oo\Ht/LiL2)Ut) . 
By Gronwall's inequality (see for example Appendix B of Evans [11), we obtain 

C^{t) = e'^^ll"!!""-^" l'"(*'/-^i^")l''*'C(0) = e'9^ll"IU(i+|in(*/i'ii2)|)t/iii2(-(Q) ^ 
But C(0) = 0. Hence C,{t) = for all t. This implies Z{t) = for all t which implies z{x,y) = for all x,y. □ 

Next we derive the explicit solution of ([8]), for the case <j> — ip = a = 1. 

Corollary 5.2. Suppose that Li,L2 > and either /? < or L1L2 < I//?- Then the unique solution of the 
Cauchy problem (0) with (f) — ijj ~ a — 1 is 

u{x,y) = (1 - I3xy)~'^ . 

Proof. Uniqueness was proved in Lemma [O] and it is trivial to check that this solves the PDE Therefore, 
assuming that u is integrable on [0, Li] x [0, L2], we may derive D'Alembert's formula by standard calculus: 



lnu{x,y) = J —liiu{x',y)dx' + lmp{y) 



(9^ /"^ d 

■\nu{x',y')dx' dy' + / —lii4i{x')dx' + lmp{y) 



(o.x)x{o,y) dxdy Jo dx 



= 2/3 / u{x',y')dx' dy' + \n^{y) +\n(f>{x) -\na. 

JiO,x)x(0.v) 

The only issue is to check integrability, which amounts to checking i'af(x,y)e[o.Li]x[o,L2] 1 ~ P^y > 0. This holds 
if and only if /3 < or L1L2 < 1/(3. □ 



Let us briefly explain one approach to deriving this formula. For nonlinear PDE's one always first guesses a 
scaling solution, in hopes of finding an explicit formula. Because of the hyperbolic nature it makes sense to look 
for a solution u{x,y) — U{xy) for some U{z). This leads to the ODE 

— lnUiz) + z—^ 
dz az'^ 

which can also be expressed as 



— lnC/(z) + z— lnC/(z) = 2pU{z), 



^ (z±\nU{z)] = 2PU{Z). 
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The idea of using a power law solution is natural because the derivative of the logarithm results in a power law, 
itself. Trying U{z) = (1 + cz)p leads to 

In [/(z) = pln(l + cz) ^ z^\n(b{z) = = p--^ => ^ ( z -^lnU{z)] = 



dz 1 + cz 1 + cz dz \ dz J (1 + cz)^ 

So, taking p = —2 and c = — /?, this solves the equation, and gives U{z) = (1 — /3z)^^ =^ u{x, y) = {1 — /3xy)^'^. 
Finally, we are led to the solution of the general Cauchy problem. 

Corollary 5.3. Suppose that i/j : [0,1] — > M are continuous and satisfy c < (f),tp < C , for some constants 
< c < C < 1. Also suppose that 0(0) — '0(0) = o; for some a. Then the Cauchy problem has a solution if 
and only if (3 < or 4>{x) dx tp{y) dy < a/p. In case a solution exists, it is unique and equals 

(a\ ( \ _ Q0(a:)-0(y) 

^ ' "^""'^^ (a - /5$(a;)*(y))2 ' 

where $(x) = (f>{x') dx' and "^{y) = ^{y')dy'. 

Proof. Suppose that u is any solution of ([H]). Let v be given by 

u($-i(ai/2a;),*-i(ai/2?;)) 



v{x,y) 



a$'($-i(ai/2a;))^"(^'-i(ai/2y)) ' 

Then, by the symmetry ([S]), f is a solution of Liouville's PDE © on the domain [0,a-^^^^{l)] x [0, ^-^/^^'(l)]. 
But v{x, 0) = ^(O, y) = 1 because $' = and 5*' = So uniqueness, the conditions for existence, and the 
formula for the solution all follow from Lemma [5TT] and Corollarv l5.2l □ 



6. Solving the Euler-Lagrange Equation 

By general principles, we know that a solution of ^ always exists: specifically, the optimizer in Proposition 
3.11 Next we calculate it, and prove uniqueness. 



Lemma 6.1. // / = g = 1, then the unique solution of ^ is given by @) for 

0(z) = 0(z) = $(z)^*(z) " ^ 



1 - 6-/9 ' ' ' ' ' \-e-P 1 - 6-/3 • 

Proof. Suppose u solves Note that lima;^o h[{x, y), (a;', y')) = d{y — y') for all (cc', y') G X. By the dominated 
convergence theorem, this implies 

lim / h{{x,y),{x\y'))u{x',y')dx' dy' ^ f f e{y - y')u{x' ,y') dx' dy' ^ f e{y - y') dy' ^ y , 
^^oJx Jo Jo Jo 



where we used the fact that u{x' , y') dx' = 1 for all y' . So 



■0(y) = limu(a;,y) = '^^ . 

Similar arguments lead to 0(x) = lim^^o "(a;, y) = -^e^^^. Since u{x,y) dy — 1 for all x e [0,1], it 
again follows from the dominated converge theorem, taking the limit a; — > 0, that ^{y) dy must also be 1. So 
Z — {1~ e"^) / (3. Checking, the reader will easily see that this gives the stated value for 0, if) and a. Integrating, 
it also leads to $ and ^. 

Uniqueness follows from uniqueness of the Cauchy problem. Corollary 15.31 Since this is the only possible 
solution, and since a solution exists, this must be it. □ 

Substituting in, and simplifying leads to the formula 

t,^, , . (/3/2)sinh(/3/2) 

(10) u(x,y) = 77. 

(e'3/4 cosh(/3[a; - y]/2) - e-f^l^ cosh(/3[a; + y - l]/2)) 
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Therefore, we arrive at the final formula. 

Corollary 6.2. As long as c < f, g < C for some < c < C < 1, and fix) dx = g{y) dy = I, the unique 
solution of ^ is 

, . ^ (/3/2)sinh(/3/2)/(a:)g(y) 

" ^' ^ cos\i{l3[F{x) - G(y)]/2) - e-P/^ cosh(/3[F(x) + G{y) - l]/2)f ' 

where F{x) = f{x')dx' and G{y) = g{y')dy'. 

Proof. Suppose that m is a solution of (|4]) under the conditions stated. Define 

"^^'^^ /(F-i(a;))ff(G-i(y))' 
analogously to the proof of Corollary 15.31 Note that F and G are continuously, strictly increasing bijections of 
[0, 1]. Using (m, we see that 

\nv{x,y) = lnu{F-\x),G-\y)) - In f{F-\x))- In g{G'\y)) 

= - \nZ-l3 f h{{F-\x),G'Hy)),{x',y'))u{x',y')dx' dy' . 

Making the change-of-variables x" — F{x') and y" = F{y'), we see that dx' — dx" / f{F~^{x")) and dy' = 
dy"/g{G-\y")). So we have 



lnt;(a;,y) = -lnZ-/3 / hiiF-\x),G-\y)),{F-'ix"),G-'{y")))v{x" ,y")dx" dy" . 
Jx 

But the Heaviside function satisfies 9(F{x) — F{x')) — 6{x — x') for any continuous, strictly increasing function 
F. For this reason, 

h{{F-\x),G-\y)),{F-\x"),G-'{y"))) = h{{x,y),{x" J')) . 
In other words, v also solves (jj]), except that 

■v{x,y)dy = / v{x,y)dx = 1, 
Jo 

using the change-of-variables formula, again. So uniqueness and the formula follows from Lemma |6. II □ 



7. Proof of Main Result 

We now explain the minor details needed to go from Proposition 13 . II to a proof of Theorem ll.il According to 
Pannes, Spohn and Verbeure's result, /i^ must converge weakly to a mixture of i.i.d., product measures, each of 
whose 1-particle marginal optimizes Si{fi, /io) ~ ^/^^C*)- But /i^ has marginals on {xi, . . . , xn) and (j/i, . . . , y^) 
equal to the product measures of f{x) dx and g{y) dy, respectively. Therefore, according to the weak law of large 
numbers (WLLN), we know that all the /l('s in the support of the directing measure for the limit of n^, must 
have X marginal equal to f{x) dx and y maginal g{y) dy. Hence, this constraint can be imposed when looking 
for an optimizer. This is actually a relevant comment because all optimizers, for all choices of a priori measure 
fiQ, have the same value/pressure: that due to the Mallows measure on Sn- For concreteness, we will now take 
/ = .9 = 1. 

Now suppose that fi optimizes the Gibbs formula. It must be absolutely continuous with respect to /io in 
order to not have the relative entropy equal to — oo. So we can write 

dt^{x,y) = u(x,y)dxdy, 
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where u(x,y) is absolutely continuous. Choosing any continuous function (j) : [0, 1] x [0, 1] — > M, with 

u{x,y)(j){x,y)dxdy = 0, 

X 

we can take 

Ue{x,y) = (l + e<j>{x,y))u{x,y). 
For |e| < l/||0||oo, we have that is a probability measure. It is easy to see that 



Si{fie,m) = '5'i(Ai,Mo) - e J (t)u\nu- J[l + e(j)]ln[l 



Since 1 + ecj) is bounded away from (and infinity) for the e we are considering, it is clear that both integrals 
above are well-defined. Moreover, it is clear that 



[1 + ecj}] ln[l + e4>\u = o(e) , 
because ln[l + e<f)] = e(/) + o(e) and J (j)u — 0. A similar calculation also shows that 

= ^,\h) + e^i'{[cf>{x,y)+cf>{x',y')]h)+0{e'). 
Since is supposed to be the optimizer, the terms linear in e must vanish: 

/ 0ulnw = - fi^{[(f>{x,y) + (f>{x',y')]h) . 
Since h is symmetric, by varying over all cj) orthogonal to u, we deduce that 

u{x,y) In u{x,y) = (3u{x,y) / h{{x,y), {x' ,y'))u{x,y) dx' dy' + Cu{x,y) . 



J X 

for some constant C . (The reason we cannot assume C = is because we left out one direction for 0, namely the 
direction parallel to u, so that there is an indeterminacy in this direction, as seen using the Riesz representation 
theorem.) In other words, we have just deduced equation ([3]). On the other hand, we have also proved that this 
equation has a unique solution given by (jlOp . Therefore, /i^ does converge weakly to the i.i.d., product measure 
of /i, where d^(x, y) = u{x, y) dx dy. 

Because of all this, if we take the empirical measure with respect to /i^, 

N 



TV 

i=l 



then this does satisfy just the type of convergence claimed in Theorem 11.11 But, taking the order statistics 
Xi < ■ ■ ■ < Xm and Yi < • • • < Y^r, we do have [xi, yi) — {X^- , 1^.) for some permutations a,T £ S^. Moreover 
(by commutativity of addition) 

N N 
i=l 1=1 

where tt = ra^^. As noted before, {Xi, . . . , X^) and (Yi, . . . , Y^) are distributed as the order statistics coming 
from Lebesgue measure, the effect of the Hamiltonian is only present in the Mallow model p^P'^''/^^ 
measure of tt. By the WLLN for the order statistics, we see that, defining 

N 

9N{x,y) = ^f{Xi,Yj)l(^(^,^iyff^,/N]{x)l{{j^l)/N,j/N]{y), 

*j'=l 
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we have that the random function giq converges in probabihty to /, everywhere in (0, 1] x (0, 1]. Therefore, since 

N , N 



i=l i=l 

we do deduce the theorem from the corresponding result for /i^. FinaUy note that taking exp{— — 1)) 
versus 1 — /3/n in the theorem does not matter, since the probabihty measures are continuous with respect to 
/3, and exp(-/3/(n - 1)) = 1 - + o(l))/n. 

8. Applications 

The ground state of the ZYg(s[2)-symmetric XXZ quantum spin system, and the invariant measures of the 
asymmetric exclusion process on an interval can be obtained from P^. See Koma and Nachtergaele's paper 
[TB] and Gottstein and Werner's paper [TJ for information about the XXZ model. For information about the 
blocking measures and the asymmetric exclusion process, we find it convenient to refer to Benjamini, Berger, 
Hoffman and Mossel (BBHM), [3]. The reader can easily deduce information for the XXZ model, since there is 
a perfect dictionary between these two. An excellent reference for this is Caputo's review [3]. 

An interesting perspective on the ground state of the quantum XXZ ferromagnet was discovered by Bolina, 
Contucci and Nachtergaele in [4]. They viewed the ground state of the quantum spin system as a thermal 
Boltzmann-Gibbs state for a classical model at inverse temperature f3 = ln(g~^). The state space they considered 
was the set of all up-right paths from (0, 0) to (m, n) E 1? (with m, n > 0). The Hamiltonian energy function for 
such a path is the energy under the path, and above the x-axis. Note that the Hamiltonian for the Mallows model 
also has a graphical representations as the number of "crossings" of the permutation. Using their representation, 
they explained some symmetries of the ground state of the XXZ model, and obtained estimates which were later 
useful in their follow-up paper, [5]. The two models are related, but only the Mallows model is manifestly a 
mean-field model. 

We consider the (nearest neighbor) asymmetric exclusion process on {1, ... , A^}, with hopping rate to the left 
•p and hopping rate to the right 1 — p, and q — (\ — p)/p. We no longer use p or pjv for the pressure, instead 
we use it for the hopping rate as expressed above. As BBHM explain, the invariant measure of the ASEP is 
a push-forward of P^. Given a permutation tt G Sn and a particle configuration r] = (771, . . . , t^^v) S {0, 1}^, 
let TTTy = (yyTTi, . . . jV-ttn)- (l'"', 0^^'^') = (1, . . . , 1, 0, . . . , 0) with k I's and N — k O's. Then, taking a random 
permutation tt, distributed according to P^, and letting 



V 



the law of TyC^^^"*^) is the invariant measure for the ASEP, with k particles and N — k holes. As BBHM explain, 
this is an instance of Wilson's general height function approach to tiling and shuffling 2ol§. 

The question we can answer is the non-random limiting density of 7y(''^^~'=) in the scaling limit, N —> 00, 
Pn = ^+ f3/4:N, fcjv — [yN\ . (Note that this corresponds to qn — I — /?(! + o(l))/A.) Namely, for a continuous 
function / : [0, 1] R, we have 



lim lim P „ 



l-P/N 



1 ^ /•! 



=1 







>e^ = 0, 



for all 2/ e [0, 1], where 



ry 

p{x;y) = / u{x,y')dy'. 
Jo 



"^Because of this, let us note that the ground state of the XXZ model is also a projection, or marginal, of the Mallows model for 
permutations (using the correspondence between the ASEP and the XXZ model [6]). This raises an interesting point for further 
consideration: are other integrable models projections of mean-field models? 
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The scaling = \ + (i/AN is the regime typicaUy called "weakly asymmetric." See, for example, Enaud and 
Derrida's paper [TD] , following the matrix method used, for example by Derrida, Lebowitz and Speer [7] . Note 
that while they considered the nonequilibrium case, we consider the particle conserving, equilibrium case. On 
the other hand, we are sure that the formula above is known. 

The integral for p(x] y) is readily evaluated. Setting 0, ■0, ^' and a as in Lemma |6. 11 

a(j){x) 



/5$(x)(a-/3$(x) *(?/')) 
^ ct){x)-^{y) 
a - I3<^{x)^{y) ' 

Substituting in, and doing minor algebraic simplifications, we obtain 

(1 - e-'^y)e-l^^ 



[l-e-P)-[l-e-f^^)[l-e-Py) ' 
From this formula it is obvious that the /3 — > limit recovers p{x] y) = y, as it should (for the symmetric case). 
Also, after further "simplifications," we obtain 

e''(5-^)/2sinh(/3y/2) 
^^^'^^ ~ e'3/4 cosh{(3[x - y]/2) - cosh{f3[x + y-l]/2)' 

In particular, one can observe that the particle- hole/reflection symmetry is manifest in this formula due to the 
invariance under the transformation (/3,x) '—^ (— /?, I — x). 

Finally, we note that we can partially undo the scaling limit by taking (3 oo with x = y + t/P (assuming 
< y < 1). Approximating smh{(3y/2) « ^e'^^/^ and noting that e"''/^ cosh(/3[a; + y — l]/2) since 
|a; + y — 1| < 1, we obtain 

p{x;y) 



1 + e* 

This is not correctly normalized due to the fact that dx = dt/j3, and (3 — > oo. On the other hand, this does 
recover the actual lattice scaling limit for the density (modulo a reflection), as has been previously calculated 
for the XXZ model by Dijkgraaf, Orlando and Reffert in Appendix A of [9]. 
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